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Displacive deformation via dislocation slip and deformation twinning usually plays a dominant role in the 
plasticity of crystalline solids at room temperature. Here we report in situ quantitative transmission electron 
microscope deformation tests of single crystal Sn samples. We found that when the sample size was reduced 
from 450 nm down to 130 nm, diffusional deformation replaces displacive plasticity as the dominant 
deformation mechanism at room temperature. At the same time, the strength-size relationship changed 
from "smaller is stronger" to "smaller is much weaker". The effective surface diffusivity calculated based on 
our experimental data matches well with that reported in literature for boundary diffusion. The observed 
change in the deformation mode arises from the sample size-dependent competition between the 
Hall-Petch-like strengthening of displacive processes and Coble diffusion softening processes. Our findings 
have important implications for the stability and reliability of nanoscale devices such as metallic nanogaps. 

Based on how atomic registries change during the plastic flow, deformation mechanisms of solids can be 
classified into three main classes: displacive, diffusional, and hybrid (coupled displacive and diffusional). 
Displacive deformations, e.g. dislocation glide and deformation twinning, tend to be prevalent at low 
temperature; the atomic registries change in a collective and semi- deterministic manner. In contrast, diffusional 
plasticity, e.g. Nabarro-Herring or Coble creep 1 , often occurs at high temperature; concurrent with the deforma- 
tion the atomic registries are randomized in an atom-by-atom manner 2 . Aside from the temperature, size scale 
can also be a decisive factor on the deformation mechanism 3 " 9 . Recently, Strachan et al. 10 found that two single- 
crystalline Au nano- electrodes (tip radius of curvature R ~ 10 nm) receded by about 10 nm after sitting at room 
temperature for 4 months. This suggests that nano -structured devices may lose their long-term shape stability due 
to diffusion as the characteristic size goes down to nanometric scale. Therefore, for the reliability and stability of 
nano -structured devices, e.g. molecular electronics 11 , nanoscale break junction 1213 , it is critical to investigate the 
size-diffusion-strength relationship. However, little has been done in the literature on this topic in a quantitative 
manner. 

In this work, we take advantage of the quantitative nanomechanical test and in situ Transmission Electron 
Microscopy (TEM) to demonstrate surface-to-grain boundary (GB) diffusional creep (Coble creep 1 ) in nanoscale 
materials. We purposefully chose Sn, a low melting point (T M = 505 K) metal, for which the homologous 
temperature (T H =T/T M ) is high (T H ~0.6) even at room temperature. This facilitates the observation of the 
size-induced switch to diffusional mechanism at room temperature, without the need for heating using special- 
ized holders in TEM or months of waiting in the mechanical tests. 

Results 

We begin with relatively large samples, with their diameters ranging from 450 to 500 nm. The compression tests 
showed that their plastic deformation process is characterized by sudden strain burst and drastic load drop 
(Fig. 1). For simplicity, we will only address the sample marked with black color in details. Kikuchi diffraction 
pattern demonstrated the tested sample which has D = 450 nm is initially single crystal. Fig. lb and Fig. lc are still 
frames extracted from the recorded movie and correspond to the points marked with b and c in Fig. la, 
respectively. The stress increased linearly with strain up to about 300 MPa (point b in Fig. la). Following this, 
a sudden strain burst occurred. At the same time, the force sensed by the loading cell dropped to zero (point c in 
Fig. la). During the strain burst, obvious geometry change occurred (Fig. lc). Slip step, characteristic of displacive 
deformation, can be identified clearly by postmortem scanning electron microscope (SEM) examination (Fig. Id). 
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Figure 1 | Compression test of Sn pillars with diameter larger than 

450 nm. (a) Engineering stress-strain curves from the compression tests, 
with average strain rate of —5 X 10~ 3 /s. (b) and (c) are TEM bright field 
images of the pillar before and after its collapse, representing the moments 
marked in the stress-strain curve for the sample with diameter of 450 nm 
(black curve), (d) Post-mortem scanning electron microscopy (SEM) 
image of the 450 nm pillar showing the slip step. 

The flat punch velocity was set to be 20 nm/s and the corresponding 
strain rate is about 5X 10" 3 /s. The data acquisition rate used for this 
test was 200 points per second which means that the interval between 
every two data points was 0.05 second. A conservative evaluation 
suggests that the maximum strain rate during the strain burst is as 
high as — 0.9/s, another indication of displacive deformation. In 
addition, the apparent yield stress (—300 MPa, point b in Fig. la) 
prior to plastic deformation (the strain burst) is remarkably high for 
Sn, since the yield strength of bulk Sn is only about 9 — 15 MPa 14 . 
The elevated strength for submicron samples is consistent with the 
well-established tenet, i.e. smaller is stronger. Similar elastic-plastic 
behavior has been reported for other submicron -sized metals 15-18 . 
Apparently, displacive mechanism still dominates the plastic 
deformation at room temperature for Sn samples in this size regime. 

Remarkable changes in deformation behavior occurred when the 
sample size was reduced below 230 nm or so, as shown in Fig. 2a. 
Significant plastic flow set in at lower yield stress (defined by loss of 
stress-strain proportionality) or even without a well-defined yield 
point. For example, the stress -strain curve of the sample with D = 
230 nm (black curve in Fig. 2a) is quite continuous except that a mild 
load drop occurred when the engineering strain reached —10% 
(point d in Fig. 2a). Fig. 2b-d are still frames extracted from the 
recorded movie that correspond to the points marked in Fig. 2a, 
respectively. A small slip step was created on the tested sample 
(marked with black arrow heads in Fig. 2d) during the mild load 
drop. For the sample with D = 200 nm (blue curve in Fig. 2a), even 
though the load drop is similar to that observed in Fig. la, the yield 
stress and amplitude of the strain burst are obviously lower. The 
lowered, blurred yield point and the large continuous plastic flow 
observed in these samples are in sharp contrast with those observed 
in larger pillars (Fig. 1), namely obvious yield point, sudden and 
dramatic load drop and strain burst. This difference suggests a 
change in the deformation mode due to the decreased sample size. 

It becomes difficult to directly FIB -fabricate samples with D < 
150 nm, because of the relatively low melting point of Sn. However, 
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Figure 2 | Compression test of Sn pillars with diameter smaller than 

230 nm. (a) Engineering stress-strain curves from the compression tests, 
with average strain rate of —5 X 10" 3 /s. (b)-(d) Still frames extracted from 
the recorded movie correspond to the events marked in the stress-strain 
curve for the sample with diameter of 230 nm (black curve). 

we found that after first compressing a larger pillar sample to "pan- 
cake" shape, the heavily deformed Sn adhered to the flat diamond 
punch. Subsequent pulling (switching from compressive to tensile 
loading) forms a narrower-D ligament, providing a means to observe 
the mechanical behavior of smaller-D samples. Compared with 
FIBed samples, these samples had fresh surfaces (no oxides generated 
in the high vacuum inside the TEM) like the mechanical break junc- 
tions 12,13 and were free of FIB damage. One of the "tensile" examples 
is shown in Fig. 3: after heavy compression, good mechanical contact 
forms between the sample and diamond tip. By pulling back the 
diamond tip, a ligament with a neck was created. Movie and images 
were taken under dark field condition using diffraction spot indi- 
cated by the white arrow in Fig. 3b. Grain boundaries (GBs) were 
observed to form through dynamic recrystallization during the afore- 
mentioned pre- deformation process: a region with bright contrast, 
meaning Bragg's law is satisfied, appeared in the neck, where 
the single-crystal ligament height (the grain size) is seen to be on 
the same order as its diameter D (Fig. 3d). As the tip withdraw, the 
necked part got thinner and thinner (Fig. 3d-e), and broke apart at t 
= 2.6 s (Fig. 3f). Interestingly, instead of forming chisel-like mor- 
phology as frequently observed in displacive deformation 19 " 21 , the 
fractured tips shrunk and receded quickly after broken. 

Subsequent quantitative tensile test of a nanoscale ligament also 
suggested that diffusional plasticity plays a dominant role at room 
temperature, as shown in Fig. 4. The probe velocity and the video 
recording rate for this test were 10 nm per second and 5 frames per 
second, respectively. Fig. 4a-d are still frames extracted from the 
recorded movie (see Supplementary Movie 1) showing the evolution 
of the sample. The starting minimal diameter of this sample (D min ) is 
about 130 nm (Fig. 4a), which decreased continuously (Fig. 4a-c) in 
response to the externally applied tensile force. At t = 5.0 second, the 
sample was seen to be torn into two parts (Fig. 4d). Again, instead of 
leaving behind chisel like sharp tips, the two parts exhibited quite 
rounded and smooth geometry. This means that in less than 0.2 
second (see Fig. 4c and 4d), the fractured sample receded into the 
geometry as shown in Fig. 4d through diffusional deformation. The 
morphology evolution is very similar to that of liquid drop. However, 
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Figure 3 | Dynamic recrystallization and grain boundary plating out 
during the tensile deformation of a nanoligament. (a) The sample profile 
prior to the tensile test, (b) Selected area diffraction pattern showing the 
sample orientation and imaging condition, (c-d) Dynamic evolution 
evidenced the formation of the new grain and deformation process of the 
nanoligament. (e) Final morphology of the fractured nanoligament. 

both the mechanical data (Fig. 4e) and high resolution TEM obser- 
vation of the sample surface (See Supplementary Movie 2) suggested 
that the samples were still crystalline solids. 

One advantage of in situ test is that the projected diameter of the 
samples can be recorded and then measured with high accuracy 
(Fig. 4e, half- filled blue circles). By assuming a circular cross section 
with minimal diameter D min , the corresponding true stress vs. time 
curve is plotted in Fig. 4e (black dots). The time points that corre- 
spond to the frames shown in Fig. 4a- d are marked with black arrows 
in Fig. 4e. It is worth noting that the apparent yield strength 
(—60 MPa) of this 130 nm sample, as defined by loss of stress-strain 
proportionality, is 5 times lower than that observed in the 450 nm 
sample (Fig. la). This "smaller is weaker" observation is very differ- 
ent from the "smaller is stronger" trend in all previous reports for 
metals (with higher melting temperatures) with similar sizes 16 ' 22 " 26 , 
and is presumably due to the significant contribution from the diffu- 
sional plasticity: atom transport along the free surface, and plating 
out at the GBs 27 " 29 . 

The strain hardening trend in Fig. 4e after t = 2.5 second can be 
rationalized by the increasing separation from atom source (surface) 
to atom sink (GB) and by the increase of the local strain rate. Along 
with the tensile deformation, necking of the sample will inevitably 
create cone-like geometry. This will generate stress gradient and 
therefore lead to localized plastic deformation. As mentioned earlier 
in this paper, the tensile test was carried out under displacement rate 
control, i.e. under fixed probe velocity. Consequently, the local strain 
rate of the necking part will keep increasing rather than being con- 
stant along with plastic deformation. This was supported by the 
increasing shrinking rate of the minimum diameter D min of the tested 
sample (Fig. 4a-c and e). By taking the derivative of D min with respect 
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Figure 4 | Tensile test of a 130 nm Sn nanoligament. (a-d) Still frames 
showing the deformation process of the Sn nanoligment. (e) The 
maximum stress plotted as a function of time (black). Instantaneous 
minimum diameters of the sample measured from the still frames were also 
plotted (blue color). Note the negative stress is due to the switch from 
compression to tension. 

to time t, we found that dD min /dt was about 5 nm/s at the beginning 
of this tensile test. However, this value increased almost 5 times just 
prior to the fracture of this tensile sample. Therefore, the strain 
hardening trend agrees well with the strain-rate sensitive nature of 
diffusional deformation mechanism. Several major stress serrations 
occurred after 2.5 seconds. This suggested that some displacive plas- 
ticity might also accompany the diffusion -dominant deformation; 
the increasing serration amplitude in true stress may indicate the 
"quantum" nature of displacive plasticity, namely, a single disloca- 
tion slip -out can cause larger and larger strain and stress relaxation 
when averaged over smaller neck volume oc D min 3 . However, no 
surface slip steps were observed during the deformation of the 
nano -ligaments. This is in sharp contrast with those observed in 
typical displacive deformation. For example, strained sub-10-nm- 
sized gold crystals always have terraces or staircases on the surface 
of the material 19 " 21 . Assuming that displacive deformation (e.g. dis- 
location burst) did take place in our test and generate terraces on its 
surface, because of the overwhelming diffusional deformation for the 
given sample size, the terraces resulted from displacive deformation 
would be smoothed out quickly by surface diffusion. 

Discussion 

Ohring et al. reported the temperature -dependent diffusivity of Sn. 
At 300 K, bulk diffusivity 30 and grain boundary 31 diffusivity of Sn are 
2.86 X 10" 17 m 2 /s and 1.45 X 10" 12 m 2 /s , respectively. Compared 
with grain boundary, free surface has reduced constraints and there- 
fore should have even larger diffusivity. Consequently, the surface 
layer may serve as a highway for short-circuit mass transport. The 
contribution of the surface diffusional plasticity to the overall plas- 
ticity should increase with decreasing sample size. Based on these 
considerations, we have constructed models to explain the observed 
diffusional deformation in this work, as shown in Fig. 5. We assumed 
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that a disordered surface layer 32 exists on the single -crystalline 
sample surface and its nominal thickness 5 S is ~ 1 nm. The surface 
chemical potential (fi) controls the drift of the surface atoms. This 
surface chemical potential is a function of external applied stress, 
surface tension, and sample geometry. If no stress is applied to the 
ligament, as shown in Fig. 5a mimicking samples in both Fig. 3 and 
Fig. 4, surface atoms will flow from the thicker parts (high fi) to the 
neck (low fi). As a result the valley is filled and the neck grows thicker 
and the whole ligament becomes more uniform, see Fig. 5a-5b. 
Fig. 5c-5d describe a possible scenario occurred in the real samples 
shown in Fig. 3 and Fig. 4. The applied stress increased the elastic 
energy in the sample, in particular the necking part, such that the 
chemical potential (fi) in the neck is even higher than the other parts. 
Diffusion flow is reversed and thus the neck is getting thinner and 
thinner. The presence of grain boundaries as a result of heavy 
deformation and recrystallization contributes to the plastic deforma- 
tion by accommodating the large plastic flow through GB diffusion 
and plating out. Specifically, the chemical potential gradient drives 
diffusional flow of atoms to plate into the GBs (see Fig. 5c-5d), 
mediating the specimen elongation under the tensile stresses. 

Tensile stress is released when the neck breaks. The fractured 
surfaces of the just broken sample have the geometry of the small 
hills shown in Fig. 5e. The hills however smooth out within just one 
video frames, 0.2 second at most. This process is driven by surface 
tension only, i.e., with no external stress at play. As shown in Fig. 5e, 

yQ 

the chemical potential of point x in the surface layer is fi(x) = — — - , y 

R(x) 

is the surface energy, Q is the atomic volume of Sn, and R(x) is radius 
of curvature. Assuming atoms in the surface layer perform random 
walk, the resulting mass flux / across the surface layer is, 
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C2 knT 



given that the mobility of atom is related to diffusivity by Einstein 
relation 33 . In Equation (1) X s is the solubility of diffusing atom in 
surface liquid layer, 3 S is the nominal surface layer thickness, k B is 
Boltzmann constant, Tis temperature, and D s is surface diffusivity of 
Sn. According to the law of conservation of mass, the deposition rate: 



d t h + CldJ s = 0. 
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Where z is the timescale of the neck evolution, we have 
h 0 yWCADs 1 



k B T R 0 



Therefore, 



h 0 Ro 3 hT 
Ty£2X s <5 s 



(4) 



(5) 



(6) 



Suppose t = 0.2s, h 0 = lx 10~ 7 m,R 0 = 5x 10 _8 m,X 5 = 1, we get D s 
~8x 10 _12 m 2 /s. This value is on the same order of that reported 
previously for grain boundary diffusivity in Sn 31 ' 34,35 , which means 
that the actual diffusion ability of Sn at 300 K is able to achieve the 
observed deformation rate. Since the boundary diffusivity limits the 
amount of mass transport, diffusion deformation is only prominent 
when sample size is small enough and the strain rate is not very high, 
depending on the material and the deformation conditions. 

The plastic deformation mechanisms of crystalline solids have 
been mapped by Ashby using the so-called deformation-mechanism 
map 36 . Depending on the applied strain rate (s) and stress (a), the 
homologous temperature (T H ) and the characteristic size scale (D) of 
the component 2 ' 37 , a specific mechanism will play a dominant role to 
carry out the plastic deformation. Normally, displacive mechanism 
will dominate at lower homologous temperatures but requires higher 
stress. Diffusive mechanism, in contrast, will only be prevalent at 
high homologous temperatures, but needs lower stress. The key 
message from our work here is that the regime boundaries of the 
deformation-mechanism map indeed shift with reducing sample 
size 2 ' 37 . 

For materials deformed through displacive mechanisms, such as 
dislocation slip or deformation twinning, it has been well established 
that the yield stress follows a Hall-Petch size-strengthening relation- 
ship, i.e. a = a 0 + kD~ a , where a is typically between l A and l 23 38 . 
The diffusional deformation, e.g. Coble creep, is strongly time- 
dependent (very high strain-rate sensitivity), without a sharply 
defined threshold stress as in displacive yielding. The strain rate 
due to Coble creep 1 is expected to be 



SCoble c 



(7) 



(3) 



where K is a dimensionless constant. Equation (7) is of the same 
origin as Equation (6). For a fixed strain rate, one can see that a 
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Figure 5 | Schematic of a diffusion dominant deformation model for the nanoligament. (a)-(b) Without tensile stress curvature driven diffusion flow 
will smooth out the necked area in the ligament leaving a ligament with more uniformed cross section (b). (c)-(d) Tensile stress increased the 
chemical potential in the thinnest part to overcome the chemical potential gradient caused by curvature. The two boundaries are sinks to the atom flow. 
(e)-(f) Sample profile evolution after fracture. 
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Figure 6 | Semi-quantitative illustration of deformation mechanism map 

of Sn. For displacive deformation, Hall-Petch strengthening predicts 
"smaller is stronger" trend (red curve). Surface diffusional deformation is 
governed by Coble creep. At a fixed strain rate, stress follows a oc D 3 , 
which shows "smaller is weaker" (black curve). These two curves crossed 
and distinguished the two deformation mechanism, i.e. displacive 
plasticity, and diffusive plasticity. 

oc D 3 . Thus, "smaller is much, much weaker" for Coble creep. Along 
with the reducing sample size, the competition between "smaller is 
stronger" displacive yield and "smaller is much weaker" Coble creep 
lead to an interesting crossover, as illustrated in Fig. 6. The intersec- 
tion of these two boundaries represents the transition between the 
two deformation mechanisms, i.e. displacive deformation and dif- 
fusive deformation. The a oc D 3 boundary will shift to the right with 
lower strain rate, or vice versa. The crossover size scale, D c (T H ,s), 
will decrease with increasing s and decreasing T H . In polycrystalline 
(nanocrystalline) metals, a transition from Hall-Petch strengthening 
to an inverse Hall-Petch relationship was predicted at very small 
grain sizes. But it has been difficult to experimentally confirm 
this expectation, because very small nanograins are unstable 
against grain growth such that samples in this grain size regime are 
usually unavailable for reliable mechanical tests. Here, using a Sn 
single crystal, we are able to demonstrate that indeed "smaller is 
much weaker", when the sample/grain size scale eventually drops 
into the regime where displacive plasticity is taken over by diffusive 
plasticity. 

In summary, by quantitative testing of Sn samples inside a TEM 
with a unique load transducer, we have demonstrated that with 
decreasing sample size the normal "smaller is stronger" trend even- 
tually switches to "smaller is much weaker". This transition arises 
from the competition between the displacive mechanisms that lead 
to Hall-Petch-type strengthening and the diffusive processes that 
mediate Coble diffusion softening (inverse Hall-Petch dependence 
on size). The dominance of diffusive plasticity as the underlying 
deformation mechanism at very small sample size is supported not 
only by the in situ observations but also by the match between the 
diffusivity calculated from an analytical model and that reported in 
literature for boundary diffusion. These findings may pave way for 
nanofabrication, for example precise sculpting of nanoscale metallic 
electrodes by surface diffusional creep at elevated temperatures, and 
have implications for life estimation of the shape stability of nano- 
metallic components at room temperature. 

Methods 

The purity of the Sn bar we used to fabricate submicron sized samples is 99.95%. For 
the compression tests, pillar-like samples with their diameters ranging from 170 nm 



to 500 nm and an aspect ratio of —2.5 were machined using focused ion beam (FIB). 
The samples were tested inside a JEOL 21 OOF TEM at 200 kV with a Hysitron PI95 
TEM Picolndenter under constant displacement rate control. The engineering stress 
a was calculated as a = 4F/nd 2 , where F is the measured force and d is the contact 
diameter. The engineering strain e was calculated as e = AL/L 0 , where AL is the 
compression displacement and L 0 is the initial pillar height. For more details see Ref. 
23. 

Metals illuminated by high-energy electron beam may suffer from radiation 
damage, including beam heating, ionization damage (radiolysis) and displacement 
damage 39 . Under normal TEM imaging conditions, the beam-heating effects are 
negligible because of the good thermal conductivity of the Sn sample (66.6 Wm _1 K _1 
at 300 K 40 ) and good thermal conduction in our experimental setup 41 " 43 . Given the 
electronic structure of metals, ionization damage is also small. Displacement damage 
is determined mainly by the threshold displacement energy (TDE) and atomic weight 
of the metal. The TDE of Sn is relatively low, at ~ 12 eV 44 . However, since the atomic 
number of Sn is large, the kinetic energy transferred from the 200 kV (the operation 
voltage of the TEM) electrons is —4.4 eV 45,46 , insufficient for displacement damage 
inside the bulk lattice. Therefore, it is the sample size induced transition from dis- 
placive-mechanism dominated plasticity to diffusional-mechanism dominated 
plasticity, rather than the electron beam effects, that plays the decisive role in the 
aforementioned change in deformation responses. 
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